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Abstract. We prove blow-up in finite time for radially symmetric solutions 
to the pseudo-relativistic Hartree-Fock equation with negative energy. The 
non-linear Hartree-Fock equation is commonly used in the physics literature 
to describe the dynamics of white dwarfs. We extend thereby recent results 
by Frohlich and Lenzmann, who established in [3| H] blow-up for solutions to 
the pseudo-relativistic Hartree equation. As key ingredient for handling the 
exchange term we use the conservation of the expectation of the square of the 
angular momentum operator. 



According to the theory of Chandrasekhar [T], white dwarfs can be described 
by a model of electrically neutral atoms interacting through classical Newtonian 
gravitation. Atoms consist of nuclei, which are responsible for the main part of 
the potential energy, and electrons, which, on the other hand, give the leading 
contribution to the kinetic energy of the star. Because of local charge neutrality, we 
can assume that the space and momentum distributions of the nuclei coincide with 
the ones of the electrons; in this approximation, we only keep track of the electronic 
degrees of freedom. Considering a relativistic dispersion E(j>) = y/p 2 + m 2 for 
electrons with mass m, and assuming a single species of nuclei with mass mz 2> m 
and charge Ze (where — e denotes the charge of the electron), this simplified model 
is described, on the microscopic level, by the quantum mechanical Hamiltonian 



where N is the number of electrons and K = Gm 2 z /Z 2 (G denotes here the gravi- 
tational constant). We use units with h — c — 1. 

Since the electron spin does not play an important role, we neglect it, and de- 
scribe electrons by wave functions in i 2 (M 3 ). In accordance with Pauli's princi- 
ple, the Hamiltonian acts then on the antisymmetric tensor product space 



In 0[6] Lieb and Yau derived Chandrasekhar equation for the ground state of the 
Hamiltonian (II. Ill in the limit of large N with nN 2 / 3 kept fixed. Additionally they 
reproduced Chandrasekhar critical number of particles N c ~ 0(k -3 / 2 ) proving 
the instability of (jl.ip for all N > N c (the white dwarf is supposed to undergo 
gravitational collapse for N > N c ). Up to a factor 4 the correct value of N c had 
already been established by Lieb and Thirring in [7] . See [3 [6j [TJ [4] for a more 
thorough discussion on white dwarfs. In principle (|l.ip can also be used to describe 
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(1.1) 




H N = f\ N L 2 {R\dx). 
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neutron stars; however, in this case, a correct understanding of the collapse requires 
the inclusion of general relativity effects. 

In the physical relevant regime of very small k and very large N, one expects 
the ground state of (|1.1[) to be approximated by a Slater determinant 

(ipi Aip2 A • •• A V>aO (x) = —t= ^2 ^V^i (2^1)^2(2^2) ■ • ■ ^n(x v n) 

of N orthonormal one-particle wave functions ipj S L 2 (R 3 ) (here the sum runs over 
all permutations of the N particles; moreover, a % = 1 if the permutation 7r is even 
while ov = — 1 if it is odd). It is simple to verify that the energy of the Slater 
determinant AjL^ipj is given by the so called Hartree-Fock functional 

(1.2) 

N 

*HF({^}jLi) = E / dx|(-A + m 2 ) 1/2 ^,(x)| 2 

3=1 J 

« f, , \Mx)\ 2 \^v)\ a -M*)4i(v)Mv)4&) 



E J dxdy 



2 .4^ J \x-y\ 

Within the range of its applicability, one also expects the Hartree-Fock theory 
to describe the time-evolution of Slater determinants. In other words, one expects 
that, in an appropriate sense, and in a suitable limit of large N and small k, 

e- lHNt (ipi A ■ • • A iPn) — Mt) A Mt) A • ■ • A ^J N {t) 

where the wave functions evolve according to the time-dependent Hartree-Fock 
equation 

N f X \ N / 1 — \ 

(1.3) idt^j = V-A + m 2 - «E [J] * ^ + * Wi) ^ ' 

Note that this system of non-linear equations, which can be formally obtained com- 
puting the variation of (jl.2[) . preserves the orthonormality relations (ipi(t), i/jj (i)) = 
5ij and the energy £hf- 

The Hartree-Fock theory can be formulated in a more compact form in terms 
of the orthogonal projection Q = Ylf=i IV'jKV'.j'I onto the subspace spanned by the 
wave functions {ipj}jLi- The Hartree-Fock energy (|1.2j) is given, in terms of Q and 
its kernel Q(x,y), by 

(1.4) £MQ) = tr vATAT^Q - J f dxdy . 

Also the time-dependent Hartree-Fock system (|1.3[) can be translated into an evo- 
lution equation for the time dependent density Q t — X^=i IV 7 .? WKV'jWI- ^ is easy 
to obtain the nonlinear Hartree-Fock equation 



(-A + m 2 ) 1 / 2 - « ( |-7 * PQ* ) + «^Q t , Q 



(1.5) ifltQ* 

where f>Q t {x) — Qt(x,x) (we denote by Qt(x,y) the kernel of the projection Q), 
and where the operator RQ t is defined by its kernel RQ t {x,y) = Qt[x,y)J\x — 
y\. By construction, it is clear that 1)1. 5|) preserves the trace TV = tr Qt, and the 



FINITE TIME BLOWUP FOR THE HARTREE-FOCK EQUATION 



3 



energy (11. 4| . In the present paper we are interested in solutions to the nonlinear 
Hartree-Fock equation (|1.5|k in particular we prove the existence of solutions to 
(|1.5p which exhibit blow up in finite time. Within the framework of the Hartree- 
Fock approximation, the blow up of solutions to (|1.5[) is interpreted as evidence for 
the dynamical collapse of white dwarfs. 

The last contribution in the commutator on the r.h.s. of (| 1 . 5[) (the term contain- 
ing the operator Rq) is known as the exchange term (while the second contribution, 
containing the density pq, is known as the direct term). The presence of the ex- 
change term is a consequence of the Pauli principle. Since, in the relevant limit of 
large N and small k, the exchange term is expected to be of smaller order com- 
pared with the direct term, it is often neglected in the physics literature. In this 
approximation one obtains the Hartree equation 



(1.6) id t Qt 



For bosonic systems (boson stars) this equation has been in fact rigorously derived 
from many body quantum dynamics in [2]. 

Recently, blow-up in finite time has been proven to occur for solutions to the 
Hartree equation (|1.6[) by Frohlich and Lenzmann in 3, 4 . To obtain this result, 
they consider the non-negative observable 

3 

M = x\] — A + m 2 x — Xj y— A + m 2 Xj 



and they estimate the expectation value tr (MQt) where Qt is a solution to (|1.6|) with 
spherical symmetry (in the sense that Q t (Rx, Ry) — Q t (x,y) for all R £ 50(3)). 
Under this assumption, they show that 

(1.7) tr MQ t < 2i 2 £ H art rC c(Q) + 0(t) 

where O(i) denotes error terms growing at most linearly in t. For initial data with 
negative energy £Hartrcc(Q) < 0, Eq. ()1.7|) leads to a contradiction to the non- 
negativity of the observable M (choosing t sufficiently large) . This implies that the 
solution Q t cannot exist globally in time. 

The spherical symmetry of the density Q plays a very important role in the 
analysis developed by Frohlich and Lenzmann; it allows them to control error terms 
arising from the commutator of M with the interaction * PQ t ) (the time 

derivative of tr MQ t contains the term tr [M, * pQ t )]Qt)- The same approach 
can be applied to solutions of the Hartree-Fock equation (|1.5|) ; it turns out, however, 
that the error terms arising from the commutator of M with the exchange term 
Rq cannot be handled like the errors arising from the direct term. This is the 
reason why the approach of Frohlich and Lenzmann does not extend in a simple 
way to the Hartree-Fock equation (|1.5[) (the method does extend to the Hartree- 
Fock equation if one assumes not only that Q is spherically symmetric, but also 
that each orbital in its decomposition Q = Ylj—i IV'j'XV'jl is spherically symmetric; 
however, as pointed out in [4], this is a physically unnatural condition). 

In our analysis, we use the strategy of Frohlich and Lenzmann; we study the 
evolution of the observable M = x \/—A + m 2 x on spherically symmetric solutions 
Qt of (|1.5p and we obtain a bound like () 1 . 7|) . The novelty of our approach lies in the 
estimate of the error terms arising from the commutator of M with the exchange 
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term Rq; to control this error terms, we make use of the expectation of the square 
of the angular momentum operator, which is a conserved quantity due to the radial 
symmetry. 

2. The main result and its proof 



The local well-posedness of the Hartree-Fock system (|1.3|) has been established 
by Frohlich and Lenzmann in [4j Theorem 1]. The well-posedness of (|1.5|) follows 
along the same line. For s > we define the space 

H s = {Q G £ 1 (L 2 (M 3 )) : \\Q\\ Hs < co} 

with the norm 

||Q||w.=tr |(1-A) S / 2 Q(1-A) S / 2 . 

It turns out that (|1 .5p is locally well-posed in TC S for all s > 1/2. 

Theorem 1 (Local Well-Posedness, [4]). Fix s > 1/2. For every orthogonal projec- 
tion Q G 7i s (L 2 (K 3 )) with N(Q) — tr Q < oo there exists a maximal existence time 
T > and a unique solution Q t G C([0, T), n s {L 2 (R 3 ))) of (Ufy with Q t > 0, and 
trQt = N for all t G [0,T). Moreover, ifT < oo, then \\Q t \\ s = tr(l - A) s Q t -> oo 
as t — * T~ (blow-up alternative) . 

For sufficiently small values of N — tr Q, Frohlich and Lenzmann also proved 
global well-posedness of the Hartree-Fock equation in [H Theorem 2] . Our goal here 
is to prove that, for sufficiently large values of N = tr Q, blow up in finite time can 
occur. In particular, we show that the time evolution of an arbitrary spherically 
symmetric density with negative energy and with finite expectation for the square 
of the angular momentum operator exhibits blow-up in finite time. To prove this 
result, we need a few simple preliminary lemmas. First of all, we need to prove 
that the spherical symmetry is preserved by the time-evolution; this follows easily 
from the local uniqueness of the solution to (|1.5I) . 

Lemma 2.1. Let Q G H 1/2 (L 2 {R 3 )), Q > 0, and assume that Q is spherically 
symmetric in the sense that 

Q(Rx,Ry) = Q(x,y) for all R G SO(3). 

For t G [0, T) denote by Q t the local in time solution to \1.5\) with Qt=o = Q- Then, 
for every t G [Q, T), we have 

(2.1) Q t (Rx > Ry) = Q t (x,y) for all R G SO (3) . 

Proof. For arbitrary R G SO(3), we define Qt{x, y) = Qt(Rx, Ry). It is then simple 
to verify that Qt is also a solution to (|1.5|) , characterized by the same initial data. 
By the local uniqueness of the solution, we immediately obtain (|2.ip . □ 

The main reason why spherical symmetry is so important to prove the blow up 
of solutions of (|1.5p is Newton's Law, as stated in the following lemma. 

Lemma 2.2 (Newton's Law). Suppose that p G L 1 (K 3 ,dx) is spherical symmetric 
with N — Jdyp(y). Then 

N 



(2-2) l«V^< N 
for a.e. x G 



A p{y) 
ay-, 



- I„.IJ 
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Proof. The proof relies on the explicit formula for radial functions 

p{y) _ i 

x-y\ \x\ 



dy ■ 



y\<\x\ 



dy ~TT- 

\ y \>\x\ \y\ 



For p e L 1 (R 3 ) n C°(R 3 ), explicit differentiation leads to 



V x / dy 



p(y) 



\x - y\ 



y\<\ x \ 



dyp(y) < t-tj. 



For general spherical symmetric p E L 1 (R 3 ) the statement follows using a simple 
density argument. □ 



An important tool in the proof of the finite time blow-up of solutions of (|1.5|) 
is the fact that the expectation of the square of the angular momentum operator 
L = x A p is preserved by the time evolution thanks to spherically symmetric 
solutions. 

Lemma 2.3. Let L = x A p denote the angular momentum operator. Let Q £ 
Ti}' 2 (L 2 (R 3 )), Q > be a spherical symmetric density with 

tr L 2 Q < oo . 

For t € [0, T), denote by Q t the local in time solution of i U.5|) with Qt=o = Q- Then 
tr L 2 Q t = tr L 2 Q for all t S [0, T). 

Proof. Observe that the angular momentum operator L generates rotations in the 
sense that 



Q V) (x) = 4>{R a x) 



for all tp £ L 2 (R 3 ). Here R a G 5*0(3) denotes the rotation around the axis d, with 
angle \a\. This implies that, for an arbitrary spherical symmetric density Q, we 
have 



ALol/ 



-iL-a 



) (x, y) = Q{R a x, R a y) = Q(x, y) 



for all a. Differentiating with respect to a, we obtain that [L, Q] = and thus that 
[L 2 , Q] — for any spherically symmetric density Q. Now, since the time evolution 
Qt of the spherical symmetric initial data Q is spherical symmetric (by Lemma l2.ip . 
it follows that 



— tr L 2 Q t = tr L 2 
dt 



VP 2 + ' 



PQ t 



nRQfQt 



tr (Vp 2 +™ 2 - k Cn * PQ?j + KRq^j [Qt,L 2 ] = . 



□ 



We are now ready to state and prove our main theorem. 

Theorem 2. (Blow-up for Hartree-Fock) Let Q € H 1 / 2 (L 2 (R 3 )) be a spherical 
symmetric orthogonal projection with N(Q) — txQ < oo, with negative energy 
£hf{Q) < 0, with finite expectation of the square of the angular momentum operator 
£ 2 {Q) = tr L 2 Q < oo, such that 

(2.3) trx 4 Q<oo, and tr (-A) Q < oo . 
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For t <E [0, T), let Qt denote the maximal local in time solution to U.5]) with Qt=o = 
Q. Then T < oo and 

\\Qtf H U2 = tr (1 - A) 1 / 2 Q t -^oo ast^T- . 

Remarks. 

• According to [4j Theorem 4] our main Theorem [2] implies that when ap- 
proaching the time of blow-up any blow-up solution exhibits a concentration 
of particles at the origin, with 

lim inf / cLt pQ t (x) > for any R > 0. 

t~*T- J\ X \<R 

• The condition £hf(Q) < requires N — trQ to be sufficiently large. 

Proof. Throughout the proof, we will use the notation p = —i\7 x . Consider the 
observable M — x\J p 2 + m 2 x — 53j=i x j\/ V 2 + rn 2 Xj. We are interested in the 
time evolution of the expectation tr(AfQ f ). 

Step 1. There exists a constant C, only depending on N(Q) and C 2 (Q), such that 



(2.4) 



— tr MQ t <tr(p-x + x-p)Q t + C. 



To prove (|2.4j) we start by computing 
d 



dt 



tr MQ t = - i tr M 



\J p 2 + m 2 — k 



PQt 



KRQ t , Qt 



(2.5) 



— i tr 
+ in tr 

— ifttr 



x\J p 2 + m 2 x, \/ p 2 + m 2 
1 



Qt 



"V 1 P 1 + rn 2 x, 



~Vp 2 + m2 x ,RQ t 



: PQt 



Qt 



Using that x = iV p , the first term on the r.h.s. of the last equation is given by 



Qt = tr (p ■ x + x ■ p) i 



— iti x\J p 2 + m 2 x, \J p 2 + m 2 
To control the second term on the r.h.s. of (|2.5p . let Vn t (x) = *PQ t ){x). Then 
i«tr 



x\Jp 2 +m 2 x,V Qt Q t 
- intY [x\j p 2 + m 2 xVq t (x) — Vq ( {x)x\J p 2 + m 2 x^j Q t 



= intr 



Ktr 



\Jp 2 + to 2 , x 2 V Qt (x) 
P 



xV Qt (x) + V Qt {x)x ■ 



\Jp 2 + 



Therefore, using Lemma |2.4[ we find 
(2-6) 

i«tr x\/v 2 +m 2 x,Vo, Qt 



< nN 



y/p 2 + m 2 ,x 2 V Qt (x) +K,N\\xV Qt (x) 



< 



K N(\\xV Qt (x)\\ + \\x 2 VV Qt \\) 
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By Newton's law (Lemma 12. 2j) . we have 



(2.7) 

and 

(2.8) 



cV Qt \\ = sup |ar| / dy 
xm 3 J 



PQ t (y) 
\x - y\ 



< CN 



;i* 2 W Qt || = sup \x\ 



dy 



PQ t (y) 



f - y\ 



< CN. 



From (|2.6p . it follows that the second term on the r.h.s. of (|2.5[) is bounded, in 
absolute value, by 



(2.9) 



in tr 



xsjp 2 + m 2 x, V Qt 



Qt 



< CkN 1 



Finally, we consider the third term on the r.h.s. of (|2.5[) . To this end, we 
decompose the density Qt in a sum over orthogonal projections 

N 

(2.io) Qt = ^\i)j,t){i)j,t\ with (i>j, t ,ipi,t) = hj- 

3=1 

It is easy to see that the wave functions {ipj,t}jLi are actually the solution of the 
Hartree-Fock system (jl.3p with initial data chosen so that, at time t = 0, 

(i/ji,ipj) — 5ij and Q = J2f=i IV'jKV'jl (it is then easy to show that the r.h.s. of 
(|2.10p is a solution to (|1.5[) ; from the local uniqueness of the solution to (|1.5|) . we 
obtain (|2.10|) ). Using this decomposition of the density Qt, we obtain 



tr jx-y/p 2 + m 2 x, Rq t 

N 



Q, 



E(V>i,i: (yjp 2 +m 2 x 2 R Qt ~ R Qt x 2 y/p 2 + to 2 ) i> ht 



3=1 



N 



(2.11) 



+ i 



xR Qt +RQ t x ■ 



N 



',;=i 



E ^ 



yV 2 + m 2 , X 2 ( r-r * l/'j.tV'i t 



3,«; 



iV 



iRe ( <A 



3,*> 



P 



1 



\ \fp 2 +m 2 

To bound the last term on the r.h.s. of the last equation, we observe that 



N 

E 


^3 




N 


< 


E 



3,t, 



VP 2 + m 2 



x ( 77 * 'fe 1 /'*,* ) V>;,t 



■ i>j,t(x) 



\J p 2 + mr 
PQtiv) 

\x - y\ 



\^i,t{x)\ \x\ 



< J dxdy PQt (x) \ X \ + E / dxd y 

< CN 2 



WM\Wi,t(v)\ 
\x - y\ 



pqAv) 
\x-y\ 



8 



CHRISTIAN HAINZL AND BENJAMIN SCHLEIN 



where we used Newton's law (Lemma 12. 2|) to perform the y-integration, and then, 
in the second term, we estimated tr(p/ yjp 2 + m 2 )Q t (p/ \J p 2 + m 2 ) < trQ t = N. 
As for the first term on the r.h.s. of (|2.11|) . we can bound its absolute value using 
Lemma HOI We find 

i,j=i L 



(2.12) 



< 



N r ( X _ \ 

,7=1 L ' 



iV 



1 



+ 



1 



To bound the first contribution, we observe that, by Lemma l2.21 



x { 77 * "<Pj,ti>i,t 



sup |a;| 



dy 



l 



r^i,i(j/)^i,t(y) 



< sup / dy 

a: J 



N - 2/| 

1 |^,t(l/)| 2 + |Vi,t(l/)| 2 



\x - 2/| 



< sup|d / dy^iM < r- y . 



(2.13) 



N - 2/1 

Next, we consider the second term on the r.h.s. of (|2.12p . For x E R 3 , we can 
find a rotation R <G 50(3) such that x — R(re 3 ) where r = \x\ and e 3 — (0,0, 1). 
Therefore 
(2.14) 



\x-y\ 6 



r I dy 



(R(re 3 ) - y ) 
\R(re 3 ) ~ 2/1 



ipj,t(y) AAv) 



3 V3 



2 

^ E 

i=l 



dy 



\re 3 - yf 



■^, t (Ry)^ itt (Ry) 



|re 3 - y| 3 

The last term on the r.h.s. of (|2.14|) can be estimated by 
(2.15) 



d VT-z — ^r? i>3-t{Ry) 4>i.t {Ry) 



<r dy 



\r~ 2/3 1 
' \re 3 - 2/| £ 



<r 2 dy 



2/3 1 



|re 3 - 2/P 



l^,t(i?2/)l 



Pt(«2/) 



< r 2 / dy 



2/3 1 



i ,3 Pt(y) 

\r - 2/1 

Introducing spherical coordinates for y — sy for y E S 2 , we find 
(2.16) 



d y — — ria V'j.i O2/) Vi.i 



|re 3 - 2/| 



< / dss 2 p t (s) / dy 



- |es-^| 



< cw. 
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Here we used the fact that, as we prove in Lemma l2~5l below. 
(2.17) 



jl-Afej 
sup / ay-. r-rr < oo . 



\>oJs 2 |e3-Ay| 

As for the first term on the r.h.s. of (|2.14p , we remark that, for example, the 
summand with i = 1 can be controlled as follows. 



(2.18) 
r 2 / dy- 



Vi 



\re 3 - y\ 



i^j,t( R y)^iA R y) 



< 



Av t. y3) f 3 ^A R y)AA R y) 



\re 3 - y\ 



r / dyy 3T 



2/1 



\re 3 - y 



ipjA R y)^*A R y) 



3^ 



Observing that 



r / dy y 3 ■ 



yi 



\re 3 - y\ 



■^jA R y) AA R y) 



r / dyy 3 d Vl 



\re 3 - y\ 



i>jA R y) ^iA R y) 



< 



r / dy {y 3 d yi - yid V3 ) 



va '>re 3 -y>^ 



+ r I dy\ yi \ } T V \ \^A R v)\\M R v)\ 
J \re 3 -y\ 6 

it follows from (|2.18p that (using the notation fj.t(y) = t/)j,t(Ry)) 

(2.19) 



r / Ay- 



yi 



\re 3 - y\ 



< 



dy ■ 



1 



+ r J dy 
< r / dy 



\re 3 - y\ 
1 



re 3 - y 
1 



ii>iA R v)^iA R y) 
(Li<Pj,t)(y)<PiAy) + 
^AmW^ A R y)\ 



r dy- 



1 



\re 3 - y\ 



(L2<Pt, t ) (y)viAy) 



\re 3 - y\ 



{\L 2 tpj, t {y)\WiAy)\ + \L2ViAy)\WjAy)\) + r / d v 



ptjy) 

\re 3 - y\ 



The last term is bounded by CN by Lemma [2721 The first term on the r.h.s. of the 
last equation, can be controlled by 



r dy 



1 



(2.20) 



\re 3 - y\ 



(\L 2 (p j:t (y)\\<p i!t (y)\ + \L 2 ^ t (y)\W]Ay)\) 



N 

< r £ / dy , (\L 2 ^, t (y)\ 2 + \^Ay)\ 2 

"HJ \re 3 - y\ 



Note that J2j=i IfjAv)? = PQt( R y) = PQt(y)- Moreover, we have 



\ L 2fjAy)\ 2 = (L 2 \ipj,t)(Vj,t\L 2 ) (y,y) = (L 2 Q t L 2 )(y,y) 
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where we defined Qt = Ylj WiAliViA- Since 

Qt{x,y) = Y^fjA^VjAv) =z2'>PjA Rx )'&jA R y} = Qt(Rx,Ry) = Q t (x,y) 



p3 



for all x,y S M 3 , it follows that Q t = Q t . Therefore, from (T2TTg]) and (|2~2T))) we 
conclude that 



r / dy-. 



yi 



\re 3 - y 



^iA R y)^iARy) 



3 ^ 



<CN + r dy 



1 



re 3 - y 



■(L 2 Q t L 2 )(y,y) 



3 

<CN + Y,r [ dy- -(L e Q t L e )(y,y). 

±1 J \re 3 -y\ 



Since ^2i = i(LgQtL£)(y,y) is invariant w.r.t. rotations of y, we can apply Lemma 
21 to obtain that 



r / dy- 



Vi 



\re 3 - y 



^iA R y)^iA R y) 



3™ 



3 

<CN + CY, dy(L e Q t L e )(y,y) 



= CN + Ctr L 2 Q t 



Using Lemma 12.31 and inserting the last bound and fj2 . 16[) back in (|2.14|) we find 
that 

- y) 



\x-y\ 



r^jAv^rAy) 



and thus, from (12.12|) . that 
I N I 

From (12.51). we obtain 



VP 2 + m 2 , a; 2 ( — * ipjjipi t 



< CN + CtiL z Q 



<CN 3 (Q) + CN 2 {Q)C 2 {Q) 



Step 2. If Qt is a solution of the Hartree-Fock equation (|1.5|) . we have 



(2.21) 

To show (|2.21l) . we compute 



— tr(p-x + x-p)Q t < 2£ H f(<3)- 



(2.22) 



— tr (p ■ x + x ■ p) Q t 

= — i tr (p ■ x + x ■ p) 



A P 2 + m 2 — K 



PQt + KR Qt > Qt 



i tr 



Now we observe that 
(2.23) 

— itr (p ■ x + x ■ p) , A P 2 + m2 Qt =2tr 



(p ■ x + x ■ p), Ap 2 + m 2 - K (n*PQtj 



nR Qt 



Qt 



Ap 2 + 



Q t < 2tr Ap 2 + m 2 Q t 
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and 
(2.24) 

in tr 



(p-X + X-p), [r-r*pQ t 



Qt = 2k tr x ■ V I i * pq 



2k I dxdyx- _ p Qt (y) p Qt (x) 
\ x V\ 

= 2k I dxdyy — — p Qt (y) p Qt (x) 



\x - y\ 3 



k / dxdy 



1 



\x - y\ 



PQ t (y) PQt 0) 



-in tr 



As for the contribution to (|2.22j) from the term with RQ t , we have 
(2.25) 

(x -p + p- x),R Qt Q t 

— 2in tr [x ■ p , R Qt ] Q t 

- 2ik tr (x ■ p R Qt Q t - Rq, x ■ pQ t ) 
'Qt(x,y)" 



2k I dx dy [ x ■ X7 X 

x-y 



= 2k J dxdy x 
= k / da; dy 



\x - y\ 
3 \Qt(x,y)\ 2 



Qt(y, x) - Q t (y, x)x ■ V x Q t (x, y) 



\x - y\ 
\Qt(x,y)\ 2 
\x-y\ 

From (222]), (|2~23]) . (|2~24|) . and (l2~25j) . we obtain (|2~2T|) . 

S'<ep 5. If Qt is a spherically symmetric solution to (|1.5p . there exists a constant 
C, only depending from N(Q) and C 2 (Q) such that 

(2.26) tr A/Q t < t 2 £ HF (Q) + < (tr (x ■ p + p ■ x)Q + C) + tr MQ . 

Eq. (|2.26p follows directly from the statements proven in Step 1 and Step 2, 
integrating twice over time. 

Step 4- Conclusion of the proof. From the assumption (|2.3|) on the initial density 
Q, it follows immediately that 



tr MQ = trxy/p 2 + m 2 x Q < tr(l + x 4 + p 2 )Q < oo 

and that 

\tr (x ■ p + p- x)Q\ < tr (x 2 + p 2 )Q < oo . 

Thus, if £hf(Q) < 0, (|2.26j) contradicts, for t large enough, the non-negativity of 
the expectation tr MQ t . This implies immediately that the maximal existence time 
T for the local solution Q t is finite. From the blow-up alternative (see Theorem [T]), 
it follows that there exists T < oo with ||Q t || w i/2 = tr (1 — A) 1 / 2 Q t — > oo as 
t-+T~. □ 

In the proof of Theorem [2l we had to control commutators of the pseudo- 
differential operator \/p 2 + m 2 with multiplication operators of the form fix) (see 
for example (|2.12p ). In this respect, it turns out that the Calderon-Zygmund theory 
of singular integrals is very useful. 
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Lemma 2.4. Suppose m > 0, p — — iV. Then, for every f € T4^ 1,00 (I 



), we 



have 



(2.27) 



vV +m 2 J(x) <C||V/|| 



A proof of this lemma can be found in [8]; see in particular the Corollary on 
page 309. The statement of this corollary does not give an effective bound on 
the norm of the commutator. However, the corollary is based on Theorem 3, on 
page 294 of [5] , whose proof provides the effective control we need (a remark in this 
sense can be found in the paragraph 3.3.5, on page 305 of [5]). 

Finally, in the next lemma we give a proof of the bound (|2.17p . 

Lemma 2.5. We have 



(2.28) 

Proof. First, we observe that 



sup 

A>o Js 



(2.29) 



sup 

A<1/2,A>2 JS : 



d y~\ — < 00 • 

|e 3 - M 



dJ x -M<c 



|e 3 - Ay| J 



because, in this regime of A, there is no singularity from the denominator. On the 
other hand, for arbitrary A e [1/2, 2], we have 



J S 2 \e 3 - Xy\ 6 



d6 smt 



1 1 — A cos ( 



2vr / d6 sin 9 



(A 2 sin 2 9 + (1 -A cos ( 
|1 - Acos6»| 

/ A 2 -l 



3/2 



2vr 



l+A 



— / dz- 

l-A 



\A 



1 — A cos ( 
2?r 



3/2 



2 A 



3/2 



\x + (l-X)\ 



3/2 



Therefore 



d y~i — - c 

S 2 \e 3 - Ay| J 



< C 



2 A 



dx ■ 



3/2 



c 



2 A 



dir. 
ZvTj2 



C 



dx ■ 



(1 + x) 



3/2 



da; ■ 



< C. 



(l-A) 



(A-l) 2 



3/2 



□ 
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